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We studied the effect of the perturbation of the meridional flow in the activity 
belts detected by local helioseismology on the development and strength of the 
surface magnetic field at the polar caps. We carried out simulations of synthetic 
solar cycles with a flux transport model, which follows the cyclic evolution of 
the surface field determined by flux emergence and advective transport by near- 
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surface flows. In each hemisphere, an axisymmetric band of latitudinal flows 
converging towards the central latitude of the activity belt was superposed onto 
the background poleward meridional flow. The overall effect of the flow pertur- 
bation is to reduce the latitude separation of the magnetic polarities of a bipolar 
magnetic region and thus diminish its contribution to the polar field. As a result, 
the polar field maximum reached around cycle activity minimum is weakened by 
the presence of the meridional flow perturbation. For a flow perturbation con- 
sistent with helioseismic observations, the polar field is reduced by about 18% 
compared to the case without inflows. If the amplitude of the flow perturbation 
depends on the cycle strength, its effect on the polar field provides a nonlinear- 
ity that could contribute to limiting the amplitude of a Babcock-Leighton type 
dynamo. 
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1. Introduction 

Surface flux transp ort models treat the evolution of the large-scale magnetic field on the 



surface of the Sun (e.g. , IWang et al.lll989t ISchrijverll200ll ; iMackay et al.ll2002t iBaumann et al. 
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2004 ). In such models, the evolution of the radial magnetic field at the solar surface is gov- 
erned by the emergence of new flux in the form of bipolar magnetic regions and by advective 
transport through large-scale flows (differential rotation, meridional circulation) and super- 
granular turbulent diffusi on. The well-known cyclic variations of diffe rential rotation in the 



form of zonal flows (e.g., iHoward fc Labontdll980l ; iHowe et al.ll2006l ) so far have not been 



considered in flu x transport simulations. On the other hand, variations in the large -scale 
meridional flow ( iKomm et al.lll993l ; iBasu fc Antiall2003l ; lHathaway fc Rightmird 120101 ) have 
been considered in flux-tr ansport simulations by assuming cycle-to-cycle changes in the over- 
all amplitude of the flow rtWang et al.ll2002a[ bikpati et alfbood ; IvVang et al.ll2009h . 



Another cycle-related modulation of the surface flow field is the modulation of the 
axisymmetric component of the meridional flow in the form of bands of latitudinal velocity 
centered on t he dominant latitudes of magnetic a ctivity, first detected at depths greater 
than 20 Mm (IChou fc Daill200ll ; iBeck et al.ll2002l ). In the case of near-surface flows, the 
residual meridional flow velocities (after subtraction of the mean flow) during cycle 23 were 
of the order of ±3—5 ms _1 and converge toward the dominant latitu des of magnetic activity 



while migrating towards the equator in parallel to the activity belts (IGizon fc Rempell 12008 



Gonzalez Hernandez et al. 2008 



motions of sunspot s and pores (Ribes &: Bonnefondlll990 



magnetic features (JKommlll994 



2010T). These flows are p robably related to the meridional 



and references therein) and other 
Meunierlll999l ). The cu mulative effect of t he near-surface 



horiz ontal flows converging towards active regions (e.g., lHaber et al.l I2004J ; iHindman et al. 



2004) appear to contribute to t he axisymmetric meridional flow perturbation (GizonJ 



2004 



Gonzalez Hernandez et al.ll2008[), but there is evidence that at lea st part of this perturbation 
is unrelated to surface activity (jGonzalez Hernandez et al.ll2010l ). 



The eff ect of the near-surface inflow s on the evolution of single active regions was recently 
studied by iDe Rosa fc Schriiyerl (120061) . Considering results obtained with a surface flux 
transport model (cf . ISchrij verll200ll ) . these authors find that inflows of the order of ~ 10 ms _1 
significantly affect the dispersal of magnetic flux from an isolated active region. These results 
indicate that the axisymmetric meridional flow perturbations associated with the activity 
belts could also affect the evolution of the solar surface field on a global scale. Particularly 
interesting in this connection is the effect on the polar field strength, which is an important 
source of the heliospheric field and also plays a significant role in Babcock-Leighton-type 
dynam o models. Here we pre sent results of solar-cycle simulations using the flux transport 
code of iBaumann et al.l ( 120041 ). including axisymmetric bands of converging latitudinal flows 
centered on the migrating activity belts. The aim of this work is to study the general effect 
of these flows on the evolution of the solar surface field, and particularly on the strength 
of the polar field. This is an exploratory study focussing on understanding the physical 
mechanisms; we do not intend to reproduce any actual solar data. We need not consider the 
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zonal flows in this study because the buildup of magnetic field to the Sun's poles is dominated 
by the latitude separation of the polarities of a bipolar magnetic region and thus essentially is 
an axisymmetric problem ( Cameron &: Schusslerll20071 ); zonal flows (and differential rotation 
in general, see ?) have no effect on the amount of signed flux reaching the poles. 



This paper is organized as follows. The flux transport model is described in Section [2j 
The relevant effects of the latitudinal flow bands on the surface flux evolution are illustrated 
with simulations of single bipolar regions in Section [3J The results of full solar-cycle simula- 
tions are presented in Section HJ which includes a study of the dependence of the polar field 
on various model parameters. The implication of our results are discussed in Section |5j 



sec 



2. Flux-transport model 

The induction equation considered in our flu x transport model is given by (for details 



Baumann et all 12004 : 1 Jiang et al.ll2009l . BOlOl ) 

dB „,, ,&B 1 d 



dt 
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dB 
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dXJR% cos 2 A d<j) 2 



(1) 



+S(X,M+D( Vr ), 

where and A are longitude and latitude, respectively, B is the radial component of the 
magnetic field, Q is the rotational velocity, v is the meridional flow velocity, rjs is the 
turbulent surface diffusivity due to the random granular and supergranular velocity field, 
S is a source term which descri bes the emergence of new flux, and the term D(j] r ) models 
the radial diffusion of the field ( jBaumann et al.l 120061 ) with the diffusivity parameter set to 
rj r = 100 km 2 s~ 1 . We use the sy nodic rotat i on ra te Q =- 13.38 — 2.30 sin 2 A — 1.62 sin 4 A (in 



degrees per day) determined by ISnodgrassI (119831 ) and take t]h = 600 km 2 s 1 



The meridional flow velocity consists of a background flow plus a perturbation, Av(X, t), 
representing axisymmetric bands of converging latitudinal flow (one per hemisphere), viz. 



v{X,t) 



v m sin(2AX) + Av(X,t) for |A| < 75° 
otherwise, 



(2) 



where v m = 11 ms 1 and 

r v sin [(A - X c {t))/AX v ] for -180° < (A - X c {t))/AX v < 180 c 
Av(X,t) = l 

I otherwise. 



(3) 
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The bands of perturbed meridional flow are characterized by their velocity amplitude, Vq, 
their width, AA„, and their central latitude, A c . The equatorward migration of the bands in 
the course of the solar cycle is represented by the time dependence of A c (see Section |4~T|) . 
Note that Equation (EJ) describes one band, its counterpart on the other hemisphere is 
obtained by changing A c — > — A c . For sufficiently small central latitudes, the two bands can 
overlap and the corresponding velocities are added. 



3. Evolution of single bipolar magnetic regions 

In order to illustrate the effect of the meridional flow perturbation on the latitudinal 
flux transport as the source of the polar field, we first study a single bipolar magnetic 
region (BMR). The temporal evolution of the corresponding surface flux depends on the 
relative position of the bands of latitudinal flow perturbation and the emergence latitude. 
We consider the evolution of a BMR that emerges at t = at a latitude of 15° on the northern 
hemisphere under the influence of four different meridional flow patterns (see Figure [T]) 



described by Eq s. (H and (H. T 



the approach of iBaumann et al.l (120041 ) 



re ini tial flux distribution of the BMR is chosen following 



Snapshots of the surface distribution of the magnetic field are shown in Figure [2J The 
four cases shown correspond to no flow perturbation (top row) and to converging flow bands 
centered on different latitudes A c . The corresponding time evolution of the polar fields 
is shown in Figure [3j When the flow perturbation is centered equatorward of the BMR 
emergence latitude (A c = ±5°, second row in Figure [2]), the overlap of the flow perturbations 
from both hemispheres (see blue curve in Figure [T]) has the consequence that preceding and 
following polarities of the BMR experience an increased latitude separation: the leading 
polarity is advected toward the equator while the following polarity is less affected. As a 
consequence, the latitudinal separation between preceding and following polarity increases, 
so that the polar field becomes stronger in comparison to the case without flow perturbation. 
The opposite effect results in the case A c = ±15° (third row in Figure [U red curves in Figs. [T] 
and|3]): now the latitudinal gradient of the meridional flow at the emergence location is such 
that the two polarities are now advected towards each other, thus reducing the azimuthally 
averaged field and, consequently, the flux reaching the pole. In the third case (A c = ±25°, 
fourth row in Figure [21 green curves in Figs. [T] and [3]), there are two opposing effects: 
the meridional flow gradient near the emergence location tends to separate the polarities 
while the following polarity experiences an overall decrease of its poleward advection, thus 
tending to reduce the azimuthally averaged field. The net effect is a slight reduction of the 
contribution to the polar field. 
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These results show that the emergence location of a BMR relative to the position of 
the bands of flow perturbation is important for its effect on the development of the polar 
field. Note that the (axisymmetric) meridional flow perturbation considered here results 
from the cumulative effect of the individual inflows. A given active region (which can appear 
anywhere in the activity belt) therefore experiences the superposition of these inflows, which 
needs not necessarily be centered on this active region. 



4. Simulation of activity cycles 

4.1. Cycle parameters 

As next step, we consider sequences of simulated activity cycles by periodically varying 
the number of BMRs that appear on the surface. The emerging BMRs have a tilt angle of 
half their emergence latitude, follow Hale's polarity rules, and are introduced in activity belts 
that migrate toward the e quator. The BMR area, A, follows the distribution n(A) oc A~ 2 



derived from observations (jSchrijver fc Harveylll994j ) . The number of BMRs emerging during 
the i th cycle is taken to vary proportional to a Gaussian time profile, viz. 

nAt) oc /exp{-[(*-^ + 6.5)/3.25] 2 } < (t - U) < 13 
\ otherwise 

where ti — 11 x i is the starting time of the i th cycle and all times are in years. With new 
cycles starting every 11 years and having a duration of 13 years we thus take into account 
the overlap of solar activity cycles. The emergence of new BMRs occurs randomly with a 
Gaussian distribution of half-width AA# about the central latitudes of the activity belts, 
A±, which migrate equatorward according to 

A±(t) = ±[Ao-(Ao-8 )(t-t,)/13], (5) 

so that the belts progress from their starting latitudes, ±Ao, to ±8° in the course of 13 
years. The resulting emergence pattern of new BMRs (butterfly diagram) for Ao = 25° and 
AA# = 6° is shown in Figure HI 

The latitudinal bands of the meridional flow perturbation move in parallel to the active 
region belts, their central latitudes (on both hemispheres), ±A C , coinciding with the centers 
of the corresponding activity belts, A±. We do not assume an overlap of the meridional flow 
perturbations from consecutive cycles; therefore, we include the flow perturbation only for 
11 years, starting from the third year of each 13- year cycle. Since the early flux emergence at 
mid latitudes affects the polar field only little, this assumption does not influence the results 
in a significant way (see also Section |4"~4"1) . 
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4.2. Dependence on the flow perturbation parameters 

Figure [5] shows the cyclic variation of the polar fields for three values of the flow pertur- 
bation amplitude: Vo = (0, 5, 10) ms _1 . The latitudinal width of the bands was kept fixed 
at AX V = 15° and the activity belt parameters were A = 25° and AX B = 6°. As already 
suggested by the results of the study of single BMRs shown in Sectional we find that the net 
effect of the flow perturbation on BMRs emerging in an extended activity belt is a reduction 
of the polar field amplitudes. The effect becomes more pronounced with increasing flow 
perturbation amplitude. 

On a more quantitative level, the dependence of the mean polar field amplitude (aver- 
ages over three consecutive cyclesjU on the width, AA„, and the amplitude, vq, of the flow 
perturbation is given in Table 1. The numbers in parentheses give the percentage change of 
the polar field with respect to the case with unperturbed meridional flow. In all cases we 
find a reduction of the polar field. For parameters roughly corresponding to the helioseismic 
results (vq = 5 ms _1 , AA„ = 15°), the flow perturbation leads to a reduction of the polar field 
amplitude by about 18% with respect to the same case but without flow perturbation. Apart 
from the reduction becoming more pronounced with increasing perturbation amplitude, it 
also is stronger for bigger AA„, i.e., for wider bands of perturbed flow. This is plausible 
because wider flow bands affect a larger proportion of the BMRs emerging in the activity 
belts and, at the same time, influence latitudinal flux advection for a longer time. In all 
cases, the evolution of the total unsigned surface flux is almost unaffected by the presence 
of the flow perturbation. 



4.3. Dependence on the activity belt parameters 

Keeping the parameters of the meridional flow perturbation fixed at values of Vq = 
5 ms _1 and A\ v = 15°, we also considered the dependence of the polar field amplitude 
on the starting latitude, Ao, and the width, AA#, of the activity belt. The results are 
summarized in Table 2. As already suggested by the results of Section [31 the biggest effect 
on the polar field occurs when BMRs always emerge near the center (latitude of convergence) 
of the bands of perturbed flow, i.e., for AX B = 0°. The flow perturbation then always tends 
to decrease the latitude extent of the BMR and thus reduces its azimuthally averaged field. 
The broader the activity belt (relative to the band of perturbed flow), the smaller is the 



1 We omit the first two simulated cycles from the analysis as these could be affected by the arbitrary initial 
magnetic field. With r/ r — 100 km 2 s _1 , the e-folding time of the magnetic field in the absence of sources 
(and thus the 'memory' of the system) is of about 5 years. 
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effect on the polar field. On the other hand, for a given activity belt width, the variation of 
the starting latitude, A , of the activity belt in a cycle does not significantly change the effect 
of the flow perturbation. In most cases, there is a tendency for the polar field amplitude to 
decrease with increasing Ao- This result may have implications for the nonlinear limitation 
of a Babcock-Leighton dynamo as further discussed in Section [5j 



4.4. Time-dependent flow perturbation amplitude 

The observations based on helioseismology indicate that the amplit ude of the axisym- 



metric flow perturbation peaks around the maximum of magnetic activity (JGonzalez Hernandez et al 



2010). We therefore also considered the effect of a temporal variation of the flow perturba- 
tion in parallel to the activity level. To this end, we modulated the perturbation amplitude, 
Vq, with the same time profile as that assumed for the number of emerging BMRs given 
by Equation (J3J), so that the maximum speed is reached at activity maximum. With the 
previously used parameters for the flow perturbation (vq = 5 ms -1 , AA^ = 15°) and for 
the width of the activity belt (AA# = 6°), this results in a polar field with an amplitude 
(three-cycle average) of 5.66 G, about 2.5% higher than the value of 5.52 G found for con- 
stant flow perturbation amplitude. The effect of the time variation is somewhat stronger if 
we assume zero spread of the activity belt (AX B = 0); in this case we obtain a polar field of 
5.22 G, which is about 6% higher than the corresponding value of 4.92 G for constant flow 
perturbation amplitude. This is to be expected since the polar field is dominated by the 
trans-equatorial transport (or cancellation) of leading-polarity flux; therefore, in the case of 
very narrow activity belts, the late phase of a cycle with fl ux emergence near the equ ator 



contributes more strongly to the strength of the polar field (ICameron Sz Schusslerl 120071 ) . 



Altogether, the effect of a temporal variation of the inflow amplitude is found to be 
rather small. Since the temporal modulation strongly reduces the flow perturbation during 
the rise and decay phases of a cycle, this result implies that the influence of the meridional 
flow perturbation on the polar field is dominated by the period around activity maximum. 



5. Discussion and conclusion 

The results presented here show that the observed cycle-related meridional flow pertur- 
bations in the form of bands migrating with the activity belts decrease the strength of the 
polar fields resulting from the latitudinal transport of surface flux. For a flow perturbation 
corresponding to the helioseismic observations, this reduction amounts to about 18% with 



respect to the case without flow perturbation. This indicates that these effects should be 
taken into account in surface flux transport simulations aiming at a quantitative pre- or 
postdiction of the polar field strength. 

It is doubtful whether this kind of flow perturbation could have significantly contributed 
to the low p olar polar field strengt h during the activity minimum between solar cycles 23 



and 24 (e.g., ISchrijver fc Liu! 120081 ) as compared to previous minima: the perturbation is 



probably present during every cycle, so that only an increase of the perturbation in cycle 23 
compared to its amplitude in previous cycles would contribute to a comparatively weaker 
polar field during the recent minimum. In any case, other effects must have been affecting 
the polar field in addition since the observed reduction by nearly a factor of 2 exceeds the 
decrease that could be caused by the flow perturbation considered here. 

The observed variation of t he flow perturbation amplitude during the activity cycle 
(jGonzalez Hernandez et al.l 120101 ) and the probable contrib ution of the near-surfac e inflows 
toward active regions to the driving of the perturbations (IGizon fc Rempell 120081 ) suggest 
that the amplitude of the flow perturbation should increase with cycle strength. According 
to our results, this would lead to a stronger reduction of the polar field built up during cycles 
of higher activity. Since, in the framework of a Babcock-Leighton dynamo, the polar field 
is a measure of the poloidal field providing the basis for the toroidal field in the subsequent 
cycle, the meridional flow perturbation is potentially important for the nonlinear modulation 
and limitation of the cycle amplitude. Furthermore, we have also seen that the polar field 
decreases for increasing starting latitude, A , of the activity b elt at the beginning of a 
cycle. Since stronger cycles typically have higher values of Ao (ISolanki et al.l 120081 ) . this 
relationship would strengthen the nonlinear effect of the flow perturbation in the subsequent 
cycle amplitude. 



W e conclude that, in addition to global variations of the meridional flow speed (IWang et al 
2002al Jbl). the cyclic perturbation of the meridional flow by converging bands migrating with 
activity belts has an appreciable effect on the build-up of the magnetic field at the polar caps. 
Its relation to the strength of a cycle means that the flow perturbation could be an important 
factor in determining the amplitude of Babcock-Leighton-type flux transport dynamos. 
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Table 1. Dependence of the polar field given in G on the parameters of the flow 

perturbation 



AA„ 



10° 15° 20 c 



■Uo[ms l ] 



0.0 


6.76 


6.76 


6.76 


2.5 


6.37 (-6%) 


6.08 (-10%) 


5.94 (-12%) 


5.0 


6.04 (-10%) 


5.52 (-18%) 


5.23 (-23%) 


10.0 


5.56 (-18%) 


4.74 (-30%) 


4.15 (-39%) 



Table 2. Dependence of the polar field given in G on the activity belt parameters 



An — 35 An — ^5 An — 15 



■u [ms l ] 





5 





5 





5 


AX B 


0° 


6.32 


4.53 (-28%) 


7.14 


4.92 (-31%) 


7.78 


5.76 (-25%) 


6° 


6.25 


5.26 (-15%) 


6.76 


5.52 (-18%) 


6.86 


5.75 (-16%) 


12° 


6.23 


5.87 (-6%) 


5.92 


5.41 (-8%) 


6.28 


5.98 (-5%) 
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Fig. 1. — Meridional flow profiles considered for the simulations of a single bipole in a time- 
independent flow. Positive flow velocities are directed northward. The black curves in all 
panels show the profile without flow perturbation (t> = 0). The coloured curves represent 
the perturbed flow profiles for Vq — 5 ms" 1 , AA„ = 15°, and different central latitudes, A c , 
of the flow perturbation. Left panel: X c = ±5° (blue curve, coinciding with the black curve 
poleward of ±20°); Middle panel: X c = ±15° (red curve, coinciding with the black curve 
poleward of ±30°); Right panel: A c = ±25° (green curve, coinciding with the black curve 
poleward of ±40° and between —10° and +10°). The dashed vertical lines indicate the mean 
latitude at which the bipolar region considered in Figure H] is initiated. 
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Fig. 2. — Evolution of the magnetic field distribution resulting from a single BMR with an 
initial unsigned flux of 1.9 10 23 Mx emerging at t — at a mean latitude of 15° and with a 
tilt angle of 5.4°. Each row shows three snapshots (after 4, 12, and 24 Carrington rotations, 
respectively). The top row corresponds to the case without meridional flow perturbation 
(vq = 0, black curve in Figure [1]). The rows below show the cases with perturbed flow 
(vq = 5 ms" 1 , AA„ = 15°) centered at different latitudes: A c = ±5° (second row, blue curve 
in Figured]), A c = ±15° (third row, red curve in Figure [1]), and A c = ±25° (bottom row, 
green curve in Figure [I]). The dashed lines indicate the equator while the full lines denote 
the central latitudes, ±A C , of the flow perturbation. Magnetic field strengths are given in G. 
The full time evolution for these cases can be viewed with aid of the animations provided 
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Fig. 3. — Temporal development of the polar cap field (average radial field poleward of ±75° 
latitude) for the simulations of single BMRs (solid lines: North, dashed lines: South) shown 
in Figure [2j The colors of the curves correspond to the four flow patterns given in Figure [TJ 
The dotted vertical lines indicate the times (^1,^2,^3) corresponding to the snapshots given 
in Figure [2J 
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Fig. 4. — Time-latitude diagram of new BMRs (+ symbols) used as input for the flux 
transport simulation of solar cycles. The red lines indicate the central latitudes, \±, of the 
activity belts. The full parts of the red lines also give the centers of the bands of meridional 
flow perturbation (beginning 2 years after the start of the corresponding activity belt). 
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Fig. 5. — Time evolution of the polar field strength (average radial field poleward of ±75° 
latitude) for different values of the perturbation amplitude of the meridional flow: v = 
(black curve), Vq = 5 ms _1 (blue curve), Vq = 10 ms _1 (red curve). The width of the bands 
of perturbed flow was taken as AA^, = 15°. 



